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GENERALIZATION OF GAS -FL0V7- INTERFEROMETRY THEORY AND INTERFEROGRAM 
EVALUATION EQUATIONS FOR ONE -DIMENSIONAL DENSITY FIELDS 
By Walton L. Howes and Donald R. Buchele 


SUMMARY 

Interferogram evaluation equations for calculating one -dimensional 
density distributions from optical-interference records are derived 
from generalized equations. The resulting evaluation equations are ap- 
plicable for any plane of focus. Assun 5 )tions involved in the deriva- 
tions are more general than heretofore and permit determination of limi- 
tations and systematic errors of the evaluation equations. Errors 
caused by an extended or misalined light source and by test-section 
windows are found to be negligible when the extension or misalinement 
is small compared with the focal length of the collimating lens. A 
criterion for applicability of the evaluation equations is established^ 
and a criterion for avoiding apparent -ray- trace crossing is derived. 

The proposed evaluation equations and procedures are compared with 
previous results and are tested in a hypothetical and an experimental 
situation involving known density distributions. 


INTRODUCTION 

Zehnder-Mach interferometers have been employed in recent years for 
quantitative studies of various aerodynamic and thermodynamic phenomena 
(refs. 1 to 6), particularly boundary layers. The optical-interference 
method possesses the in 5 )ortant advantages of permitting instantaneous 
quantitative recording of the entire phenomenon without disturbing it. 
However, because the interference method is integrative, it is presently 
restricted to certain quantitative applications by mathematical limita- 
tions that are primarily attributable to the geometry of the phenomenon 
and to optical refraction. 

Practical interferogram evaluation equations for one -dimensional 
density fields, where the density gradient is essentially perpendicular 
to the incident light path, have been proposed in references 7 to 11 
and include the effects of optical refraction. In references 7 to 10 a 
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sin5)le evaluation procedure is developed; however, this procedure in 
its most general form necessitates a "shifting" procedure, becomes in- 
valid in the presence of apparent -ray-trace crossing, and may not apply 
in the immediate vicinity of surfaces. In the present report the analy- 
sis given in reference 11 is generalized and modified in an attempt to 
eliminate the preceding difficulties and to determine the importance of 
various systematic errors. 


GENERALIZED EQUATIONS 

General functional relations giving the refractive -index distribu- 
tion as a function of measurable quantities and considering a laterally 
extended light so\irce are derived. For subsequent simplicity, the fol- 
lowing assumptions are made; 

(1) The interferometer optical system is perfect. 

(2) The light incident at the test section, which contains the 
phenomenon to be studied, is perfectly collimated. 

(3) The refractive index of the medium to be studied is very nearly 
unity. 

(a) The test medium possesses no internal space discontinui- 
ties of refractive index, such as shock waves, in the region of 
interest. 

(b) The test medium is bounded by plane discontinuities, such 
as wind-tunnel windows, perpendicular to the optical axis of the 
interferometer optical system. 

(4) All other media traversed by the light possess a constant 
refractive index. 

References 9 and 11 have shown that the pertinent effects of wind- 
tunnel windows are insignificant for the case of an axial point-light- 
source. It will presently be shown that window effects are also essen- 
tially negligible with respect to light emitted from off-axis source 
points. Hence, the effects of windows bounding the test section will 
be disregarded temporarily in the subsequent analysis. 

The inportant optical quantities involved in evaluating an unknown 
density field from an interferogram are: fringe shifts, which denote 

observed changes of phase and order of interference associated with 
density changes at fixed points in the field; and, a form of optical 
distortion (very similar to a mirage) that results from apparent shifts 
of location of points in the field. 
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In order to express fringe shifts and distortions as functions of 
measurahle quantities, the light path must he known. Light-ray traces 
in an optical medium follow Fermat ^s principle, namely, 

\|r=c y*dt=y*|jLds = extremum (l) 

which defines an optical -path length \|r (all symbols are defined in 
appendix A) , where c is the velocity of light in vacuimi, t is the 
time required for the light to traverse distances measured along the 
light paths, and \i is the refractive index of the medium. Fringe 
shifts N are defined by 

N = - ^^)/\ (2) 

where X ^denotes a specific wave length of light, and optical lengths 
and are associated with coherent light waves and are measured 

from the light source to the plane at which interference is recorded. 

Two or more light waves are coherent and, hence, may be observed to in- 
terfere at the interferometer image plane when any pair of ray traces 
associated with the waves proceeds from the same light-source point and 
intercepts the same image point. 

A few representative ray traces throughout an interferometer opti- 
cal system are shown in figure 1. Interference at a given image point 
may be represented generally in terms of one ray trace contained in the 
test beam (trace l) , which traverses the test section, in conjunction 
with a second ray trace contained in the reference beam (trace 2), which 
circumvents the test section. However, it can be easily demonstrated 
(ref. 8) that those fringe shifts which axe of practical interest may be 
expressed in terms of ray traces in the test beam alone. Thus, the 
word "reference” will henceforth refer to reference conditions pertain- 
ing to the test beam. The fringe shifts of immediate practical interest 
are : 

(1) N^, associated with the refractive index change to 

within the test section, where the subscript a refers to 
atmosphere external to the test section, and the subscript * 
refers to ambient conditions within the test section 

(2) N, associated with the refractive index change to q with- 

in the test section, where generally q = q (x,y,z) and x,y,z 
are cartesian space coordinates 

Let the x,y,z coordinate system be a right-handed system, as de- 
fined in figure 1. Then, including fringe shifts contributed by off- 
axis points in the light source, it is shown in appendix B that 
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N = : 



( 3 ) 


where, as shown In figure 2, L is the test-section span, <»> is the 
angle that a particular ray trace makes with the z-axis, D is the opti 
cal distortion which has "been defined previously in reference 11 and 
vill subsequently he discussed, and K is a fraction of the test- 
section span defined by 


where the real object-plane z = Zj^ is the real locus of objects imaged 

at the selected image plane. All geometrical quantities in equation 
(3), except those quantities bearing the subscript r, are referred to 
ray traces that traverse the field p. The subscript r refers to ref- 
erence traces that traverse the field p«>. Quantities denoted (-1) or 

(+l) are referred to the medium immediately preceding or succeeding the 
test section (atmosphere, as a result of the present assun^jtions) , re- 
spectively. Equation (3) represents a set of equations for each point 
of the interferogram, one equation corresponding to each pertinent pair 
of ray traces which intercepts a given image point. 

If Pg^ replaces p^ as a reference and p reduces to the constant 
val\je p , then equation (3) reduces to 


associated with an axial point-light source (ref. 12). 

Distortion of the locations of conjugate points can be derived by 
inspection of figure 2. Thus, 


Ng = -[p.L sec <P - p^L sec - p^D sin <p (-l)] 
X. 

which is mosh commonly applied in the paraxial form 


(3a) 


Nc = ^ L 


(3b) 



where 
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Zl = ±|(Z - + 211 ^ [l - I 

Z = KL I sec <P(+1)| 

Z^ = KL I sec <P^(+1)| 


COS V 




(5a) 


(5b) 

(5c) 


and V is the angle of refraction, H is the refraction displacement 
perpendicular to the z-axis at plane z = L, and u is the angle formed 
hy the intersections of a trace through the field [i, its unrefracted ex- 
tension from plane z = 0, and the associated reference trace through 
the field at, and in, the plane z = L. 

Eqimtions (3) and (s) are \aseful for starting the derivation of 
interferogram-evalmtion equations for particular refractive-index- 
field geometries and for investigating effects of extended or off-axis 
light sources. 


Becavise of mathematical con 5 )lications, only the following cases 
have been considered: 

(l) |i = constant 


An evaluation equation for case (l) is supplied, in general, hy equation 
(3a) and for the particular case of the axial point -light -source by 
equation (3b). Evaluation methods for case (2) have been described in 
references 7 to 11. Evaluation methods for case (3) have been described 
extensively in several papers including references 1, 13, 14, and 15. 

The present report is concerned with case (2) in which |jl is a func- 
tion of a single cartesian coordinate, say y, which is perpendicular 
to the optical axis. (The density p is related to the refractive 
index by the Dale -Gladstone formula 


where k, the specific refractivity, is a constant for any given wave 
length.) However, it is noteworthy that equations (3) and (5) also re- 
duce to the standard starting equations employed for evaluating axially 
symmetric distributions M’(r) (refs. 1, 13, 14, and 15). 


APPLICATION TO ONE -DIMENSIONAL DENSITY DISTRIBUTION 


(2) p = p(y) or p(x) 

(3) p = p(r) (r^ = x^ + y^) or 


p = p(r) (R^ = x^ + y^ + z^) 


fi = 1 + kp 
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Evaluation equations for density distributions of the form p(y) 
are derived in series form in appendixes C and D by utilizing equations 
(3) and (5). For simplicity the evaluation equations are expressed in 
terms of the object-space coordinate y, which is related to the cor- 
responding image -space coordinate y' by 

y’ = -my 

where m is the lateral magnification introduced by the camera lens. 
Also, the distribution p(y) is evaluated in terms of the density ratio 
p*, where at a given ordinate value y^ 

p(y„) 

p* . 

The present evaluation equations are based, whenever possible, upon 
sin^jlifying assumptions involving the independent variables of an ex- 
periment, for example, refractive index, geometry of the phenomenon, or 
geometry of the experimental apparatus, rather than, for example, the 
light path. The principal additional assumptions are as follows; 

(1) The refractive index p is a monotonic function of y along 
any given ray trace. 

(2) The function p is representable along any given ray trace by 

the power series p^ = where the cartesian coordTnate t) is 

defined by tj = |y - yQ|. This permits sinplif ication and extension of 

the analysis to decreasing, as well as increasing, functions p(y). 
Approximations are denoted one-term, two-term, ... according to the 

number of terms of the series which is utilized. 

(3) Ray traces are representable by the power series t] 
in the interval o = ^ = L, where ^ = z. 

According to the preceding assumptions, one- and two-term approximations 
are associated with the following additional assumptions: 

One -term approximation: 

(a) Along any given ray trace, p = constant. 

(b) The ray traces are straight lines through the test section. 

(c) The measured fringe-shift profile is always distortionless. 
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(d) The measured fringe -shift profile is a step function of 
y or a constant for all values of y. 

Two-term approximation: 

(a) Along any given ray trace p is a linear function of y. 

(h) The ray traces are parabolas within the test section. 

(c) The measured fringe-shift profile is distortionless 'when 
the selected object plane corresponds to the midspan plane of the 
test section. 

One- and two-term approximations and the corresponding evaluation 
procedures associated with an axial point-light-source are as follows: 


One -term approximation: = Iq 


(a) 

P* 

= 1 + 

N : 

Two-term 

approximatl( 

(a) 

p* 

= 1 + 

_1 

kp, 

(b) 

7 = 

^0 = 

yp 

(c) 

D = 

- - (1 ■ 
2 '' 

(<i) 

bl 

II 

X 

L 




N 


^T1 0 

X 1 


£ - f (2 - 3K) l2J 


— (incidentally, b^ = k 


(i) 


yo 


, also. 


Evalmtion procedure: 

(a) Plot N as a fimction of graphically from measured 

data. 

(h) At each datum point measure the slope of the resulting 
profile Nj). 

(c) For each dattim point compute then p*, I), and jq 

(if necessary). The density ratio p* as a function of 
y represents the desired profile. 


The convention + or - is associated vrLth p(y) increasing or de- 
creasing, respectively. 
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Although, when K = l/2, the preceding one- and two-term evaluation 
equations are identical to the corresponding results indicated in ref- 
erence 11, important differences are made evident by the derivations, 
namely: 

(1) The present evaluation equations apply for any value of K and 
permit direct calculation of p* without recourse to a shifting pro- 
cedure (refs. 8 to 10). 

(2) The additional assumption dD/dy = 0 (stated somewhat incom- 
pletely as D = 0 in ref. ll) is unnecessary (but instead follows as a 
direct result of the two-term assumption), so that mention in reference 
11 of an iteration process for eliminating residual distortion is in- 
consistent with respect to the primary analysis. 

By relating derivatives of the measured fringe-shift profile to 
derivatives of the corresponding hypothetical distortionless profile 
and then to coefficients b^, it follows from the derivations in appen- 
dixes C and D that, as a result of the two-term approximation, 
dN/dyp = dU/dy because do/dy = 0; that is, the slope of the measured 

fringe-shift profile at y = yj) has the same value as that of the cor- 
responding hypothetical distortionless -fringe-shift profile at y = Jq. 
Thus, correct values of the coefficient b^^ are presumably determinable 

from the measiu'ed fringe -shift profile associated with any value of K, 
and, hence, the present equations apply for any value of K. 


Theoretical Limitations 


The analysis described in appendixes C and D can be used to obtain 
theoretical information concerning limitations of the interferometric 
method and the preceding evaluation eq\iations and concerning the im- 
portance of certain other systematic errors. In the succeeding para- 
graphs the following items are considered: 

(1) Establishment of a criterion for minimizing apparent -ray- trace 

crossing 

(2) Series convergence and the calculation of remainders associated 

with the finite approximations 

(3) Establishment of a criterion for the applicability of the pro- 

posed evaluation eqmtions 

(4) Estimation of the importance of lateral extension or misaline- 

ment of the light so-urce 
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(5) Determination of desirable light-source geometries 

(6) Estimation of the effect of test-section windows with respect 

to off-axis light-source points 

Apparent -ray-trace crossing , - Apparent -ray-trace crossing, which 
has been noted in references 9 and 10, occurs when the relative magni- 
tudes of the initial ordinates Yq and corresponding apparent ordinates 

Y-£^ associated with any pair of refracted ray traces are reversed (as in 

the case of an inverted mirage). Two cases of apparent -ray-trace cros- 
sing are illustrated in figure 3. In figure 3(a) ^^crossing” prevails 
throughout the interval for which 



whereas in figure 3(b) crossing occurs throughout the interval for which 



It is demonstrated in appendix E that crossing occurs whenever 


Ay 


^ 1 


Regions of an interferogram that are associated with crossing ex- 
hibit multiple imagery in that records of two or more regions of the 
field p(y) are superimposed and interrelated. Such regions probably 
cannot be evaluated unless one trace of each of the crossed pairs is 
prevented from reaching the final image plane. Unfortunately, the ex- 
istence of crossing may not be apparent from the interferogram. 


The possibility of crossing can generally be alleviated by the 
proper choice of K. It is shown in appendix E that, for arbitrary 
functions p(y), apparent -ray- trace crossing is least likely to occur 
when K == l/s because do/dy is essentially zero when K = l/2. A 
criterion for avoiding apparent -ray-trace crossing is also derived in 
appendix E as follows: The existence of apparent -ray-trace crossing is 

unlikely when the inequality 


± i (1 - 2K) k ^ l2 



(1 - 4K) 




^ 1 


is satisfied. The inequality may he satisfied by properly selecting 

K if information regarding the expected extreme value of d^p/dy^ is 
available prior to the experiment. 
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Series convergence. - It has been shown In reference 9 that, if the 
series expansions for p, and t] converge, then the resulting power 
30 P 10 S for p and D also converge. The xntroduction of a laterally 
extended light source still yields a convergent series, hut differentia- 
tion of series may not always yield a convergent series. In an arbitrary 
experiment the existence of convergence is not assured, although it 
would appear that the power series for p* and D are usually con- 
vergent or, at least, asymptotic because the one-term expression for p* 
has always been found to yield results which are in reasonably good 
agreement with predicted physical values of p*. 

The remainders associated with the one- and two -term approximations 
of p* and D are given in appendix E, where remainders associated 
with the two-term approximation depend upon the three-term approximation 
which is derived in appendix F. 


Applicability criterion . - A criterion for the applicability of the 
two-term approximation is derived in appendix E. It is shown that the 
two-term approximation is likely to be valid for at least those regions 
of the field for which the second derivative of the measured fringe- 
shift profile satisfies the inequalities 


XL 


< ± 


-1.7 

XL 


< + 


- 0.8 

XL 


< 


± 


d^W 


dy 


D 


ii 

XL 


when 


d^N ^ 1^ 
dy2 XL 


when 


d^N < ^ 
dy2 XL 


when 




Calcuilation of p* and D up to and including the two-term ap- 
proximation is practical. Higher -order approximations are excessively 
conpllcated and are limited by the inaccuracy of measuring higher -order 
derivatives of N. Also, no practical situation has yet been en- 
countered in which the additional terms involved in the three-term ap- 
proximation are appreciable. 

Extended or misalined light source . - Estimation of the importance 
of off-axis light-source points resulting from lateral extension or mis- 
alinement of the light source can be made by inspecting the equations 
for p* and D which are derived in appendix C. Thus, if p and q 
represent coordinates of light-source points in a plane perpendicular 
to the optical axis such that 
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p = - X 4- constant 




where the origin (p^q) = (0,0) lies on the optical axis^ then the ef- 
fects of off-axis source points are expressed by 


where 6 denotes the variation of value of the associated quantity from 
its value when the light source consists of an axial point, and f is 
the focal length of the collimating lens. The pertinent quantities are 
indicated in figure 4. Experimental evidence, which qualitatively con- 
firms the preceding conclusion, is presented in reference 8. 

Light -source geometry . - It is found by inspecting the series ex- 
pansions for p* and D in appendix C that, from the standpoint of 
evaluation -equation sinpilicity, light -source geometries in order of de- 
creasing desirability are: 

(1) Axial point 

(2) Narrow slit along p-axis 

(3) Narrow slit along q-axis 

(4) Extension in p,q-plane 

Test-section windows . - For simplicity, the effects of test-section 
windows upon N and D were disregarded in the preceding analysis. 

It is shown in appendix E that window effects are negligible when 


CD 



and 


5D « 0 



<<1 



and 
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Analytical Verification 

The validity of proposed evaluation procedures can he investigated 
and their limitations illustrated by calcvilating exactly the profile 
Nj) associated with an ass\med hypothetical profile p(y) and then by 

attenpting to recalculate p*(y) from the proposed evaluation equations 
and the computed profile Np. However, the equations for computing Np 

m\ast be assumed correct and exact, or nearly so. Otherwise the evalua- 
tion equations may effect an excellent reproduction of the assvimed pro- 
file p**(y) although they are incorrect, as was the case in reference 16 

The present evaluation equations were tested by assuming an expo- 
nential profile p(y). The profiles Np were calc\ilated accurately and 

the recomputed profiles P*(y) were calculated to one-, two-, and three- 
term approximations for three values of K, namely K = l/3, l/2, and 
1 and for two values of L, namely L = 1.8 and 3.6 inches. The equa- 
tions for calculating Np are given in appendix G. The recomputed pro- 
files p*(y) for L = 1.8 and 3.6 inches are presented in fig-ures 5 and 
6, respectively. The two-term profiles for K = l/2 and L = 1.8 and 
3.6 inches are compared in figures 7(a) and (b), respectively, with the 
corresponding profiles obtained by using the Wachtell-DeFrate method 
(refs. 7 to lO). 

A value of the exponential coefficient a (appendix G) was selected 
that would introduce strong nonlinearities in the profile p(y) and 
hence permit illustration of limitations of the interferometric analy- 
sis and evaluation equations. 

For L = 1.8 inches, the two- and three-term approximations yielded 
recomputed profiles P*(y) that are in excellent agreement with the 
assumed profile p*(y) for all values of K (as suggested in the pre- 
vious discussion) and represent definite improvements ovei the results 
of the one -term approximation. It is apparent from fig\jre 7(a) that 
the two-term approximation for K = l/2 permits evaluation closer to 
the surface (y = O) than does the Wachtell-DeFrate method. In particu- 
lar, p(y) is only determinable for 

y ^ 0.0032 inch for K = l/3 (method of refs. 7 to 10) 
y = 0.0012 inch for K = l/2 (method of present report) 

For K = 1/2 and 1 the maximum absolute magnitudes of d%/dyp slightly 

exceed the values set by the applicability criterion, but the agreement 
of the two-term res\ilt with the assumed profile p*(y) is still very 
good. 
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The ratios of successive terms of the series expansions for p* 
and D are proportional to L. Therefore, as L is increased, higher- 
order terms increase in relative importance. For L = 3.6 inches, the 
two-term approximation represents an improvement over the one-term ap- 
proximation, and the profiles for all values of K are in agreement; 
however, because of the increased value of L the y interval of good 
agreement with the assumed profile p*(y) is considerably reduced. In 
particular, the recon5)uted profile p*(y) for K = l/3 is strongly af- 
fected by apparent -ray-trace crossing and the recomputed profile p*(y) 
for K = 1/2 is slightly affected, as is apparent from figure 8. How- 
ever, distortion is smallest when K = 1/2 and is an absolute over -all 
minimum for some plane in the interval 1/2 < K < 1, as is |dl)/dy|. 

For all values of K, the maximum absolute values of d^N/dy^ exceed 

the values set by the applicability criterion. However, the profiles 
P*(y) recomputed according to the two-term approximation do not begin 
to deviate appreciably from the assumed profile until the criterion has 
been definitely violated. Thus, it may be tentatively surmised that 
the applicability criterion is somewhat conservative. The three-term 
profile P*(y) is strongly affected by the preceding difficulties and 
is less useful than the profile obtained from the two-term calculation. 
According to figure 7(b) the Wachtell-DeFrate result is again limited 
to a smaller y interval than the present two-term result. 


Experimental Verification 

In order to verify conclusively the existence and magnitude of the 
optical-refraction effect, an experiment was devised for measuring the 
effect and confirming the refraction terms associated with the two-term 
approximation of p* and D. Good agreement was obtained between 
theory and experiment. Details and results of the experiment are dis- 
cussed in appendix H. 


Comparison with Other Analyses 

Objections to analyses presented in references 16 and 17 have been 
discussed in references 9 and 10. Objections to the analysis in refer- 
ence 11 have already been discussed in the present report. 

Results given by the two-term equation derived independently by 
Wachtell and DeFrate (refs. 7 to 10) are equivalent to those given by 
the present one -term equation when K = 1/3. The former equation may 
be readily derived from the present two-term equation. For any other 
selected value of K, Wachtell in references 9 and 10 proposed a 
** shifting" procedure by which an uncorrected density profile is com- 
puted and then corrected by means of shifting equations to give the 


14 


NACA TN 3340 


profile associated with K = l/3. However^ if apparent -ray- trace cross- 
ing occurs when K = l/3^ then the entire method fails because the den- 
sity profile must be corrected to K = l/3. For K = 1/3, 


so that a ray trace that satisfies y = 0 at z = 0, where y = 0 repre- 
sents the surface of a model, will generally appear at y = D / 0 (ref. 
8). Thus, the Wachtell-DeFrate method can only be used to evaluate 
in the immediate vicinity of a surface when (dp/dy)o = 0* 

By comparison, failure of the present method resulting from 
apparent -ray-trace crossing should not occur beca-use the corrected den- 
sity profile is associated with the selected value of K, rather than 
with K = 1/3, specifically. Generally, a value of K can be chosen 
for which crossing is absent. Moreover, the computed distortion is zero 
if K = 1/2. However, slopes dN/dy^ must be measured, and a refrac- 
tion term must be computed in calculating p*. On the other hand, 
p*(y) can be computed for any other selected value of K without re- 
course to a shifting procedure if D is computed. The total amount of 
con^jutation is equivalent to that of the Wachtell-DeFrate method with 
shifting. 

In reference 8 an iteration procedure is proposed which presumably 
can be applied when p(y) is extremely nonlinear and when there is no 
apparent -ray- trace crossing. However, the iteration method in compari- 
son with the preceding approximation methods necessitates an excessive 
amount of computation and some degree of guesswork. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 19, 1954 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report : 


a 

exponential refractive-index-distribution coefficient 

h. 

refractive-index coefficients 


by 2 refractive -index coefficient; first and second approxima 

’ ’ tions, respectively 


c 

constant 

c 

velocity of light 


light-path coefficients 

c* 

» a 

light -path coefficients 


fringe -shift coefficients 

D 

optical distortion 

e 

exponential 

f 

focal length of collimating lens 

H 

refraction displacements of ray traces in and at plane 
z = L 

K 

Zp 

fraction of test-section span, K = 1 - -rp 

Li 


constant, = t - aK 

k 

specific refractivity 

L 

test -section span 

It ,1,1 

lengths given by equations (5a), (5b), and (5c), respec- 
tively (fig. 2) 

m 

lateral magnification 

N 

interference -fringe shift 
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function representing hypothetical distortion-free 
interference -fringe -shift profile 

function representing measured (distorted) interference- 
fringe-shift profile 

integer 

light -source coordinate parallel to x-coordinate 
light -source coordinate parallel to y-coordinate 
radius in spherical coordinates 
remainder 

radius in cylindrical coordinates 
length measured along ray trace 
absolute temperature 
time 

window thickness 

quantity defined in appendix G 

angle formed by intersection of refracted-ray trace , its 
unrefracted extension from plane z = 0, and its associ- 
ated reference trace at, and in, plane z=L (fig. 2) 

refraction angle of refracted ray-trace (fig. 2) 

cartesian coordinates perpendicular to optical axis 

y-coordinate value of ray trace at plane z = 0; y-coordinate 
values associated with profile Nq 

apparent y-coordinate value of source of ray trace; 
y-coordinate values associated with profile 

axial separation of real- and apparent -object plane 

optical axis 


real-object plane 
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a 

projected angle in xy- or -plane 

P 

projected angle in yz- or -plane 

r 

projected angle in xz- or 4?-plane 

A 

increment 

6 

variation or change 


ray-trace coordinate parallel to optical axis 

n 

ray-trace coordinate perpendicual to optical axis 


coefficients 

X 

wave length of light 

P 

refractive index 

V 

integer 

4 

ray-trace coordinate perpendicular to optical axis 

p 

density 

p* 

density ratio, p* = p(yo)/p« 

a,T 

integers 

<P 

angle formed by ray trace and optical axis 

t 

optical -path length 

(-2),(-l) 

refers to mediums preceding test section 

(+l),(+2) 

refers to mediums succeeding test section 


Subscripts : 


A 

refers to lower siirface of upper block in hot-plate model 

a 

atmosphere 

c 

refers to change of refractiv-e index within test section 
from that of the external atmosphere to some other con- 
stant value 


18 


NACA TN 3340 


g window 

r reference ray trace 

w value at wall, or surface 

m,n,v,a,T integers 

“ refers to ambient conditions within test section 

Superscript : 

’ refers to image space associated with camera lens 


maximum absolute value 
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APPENDIX B 

GENERALIZED FRINGE-SHIFT EQUATION 

In general, the interference fringe shift measured at a given image 
point is represented by equation (2), namely 

N = (?i - i2)A A) 

The quantities and are determined from associated ray traces 

whose paths in any given optical medium are found from the first-order 
partial differential equation (ref. 18) 



Equation (B1) is the solution of the variation problem indicated by 
equation (l), where q = |-l (x,y,z) may be an arbitrary function of the 
cartesian coordinates y, z. 

If and ^2 associated with ray traces which traverse the 

fields |JL and respectively^ then Aijr = - ^1^2 is determined by 

ray traces 1 and 2 in figures 1 and 2; these traces satisfy the re- 
quirements for interference because they are parallel in the space pre- 
ceding the test section. It is shown in references 7, 8, and 11 that 
A\(r is determined entirely in the interval bounded by the plane wave- 
fronts in the medium (-1), which immediately preceeds the test section, 
and by the effective plane, or spherical, wavefronts determined by 
traces 1 and 2 in the medium (+l). Thus, N may be expressed by inspec- 
tion in terms of the geometrical quantities shown in figure 2. The re- 
sult is equation (3), in which the first term on the right-hand side is 
attributed to the trace through the field p, the second term to the 
trace through the field p^, the third term to obliquity of traces in 
medium (-1), and the last term to obliquity of traces in medium (+l) 
plus refraction. 
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APPENDIX C 


DERIVATION OF EVALUATION EQUATIONS ASSOCIATED WITH A 
ONE -DIMENSIONAL DISTRIBUTION p(y) 

Light Path 


Light ray traces are determined by 



For a function |jl = \i{y) ^ equation (B1) possesses a solution of the form 

\lf(x,y,z) = + ^ 2 (y) + ^ 3 ( 21 ) 

namely 


Ot 


t = ± C-i (x - Xq) ± 


yo 


- c| dy ±C 3 (z - zq) 


where 


'1 2 _ p2 


cv + 


are constants. Applying Jacobies theorem (ref. 18) reduces the solution 
for 'll; to solutions for the ray traces, namely 


dz 


= ± ^ - c| = tan p 


g = ± 1 = tan 


a 


Initially, when 


X = Xr 


y = yo 


Z = Zq = 0 


then 


a = a,(x^) 


P = P(o) 
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and the constants Cj and Cj are given by 

12 


C? = 


tan^ y(0) 


^ 1 + tan^ P(o) + tan^ t(0) 


C# = 


[r-(yo)]‘ 


1 + tan^ 3(0) + tan^ t(0) 


where 


dx 

dz 


tan X 


Thus^ the ray traces are given hy 


1 .(gf . tan^ r(0) = 


ti(y) 

>^(yo). 


sec^ ‘P(O) 


where 


sec^ 4>(0) = 1 + tan^ P(0) + tan^ t(0) 


and 


(Cl) 


^ = constant 

Equation (Cl)^ which applies when |i(y) is an increasing function^ 
can he extended to include cases where M.(y) is a decreasing function 
and, also, ultimately to lead to simpler evaluation equations by intro- 
ducing aiixiliary, floating, coordinates ? defined by 


C = |x - XqI 

n = |y - yol 

^ = z 

The coordinates of any given ray trace are then = (0^0,0) at 

z = 0. In terms of the new coordinates, |a^ corresponds to |J.(y) and 

^^0 - ^^(yo) “"O = ^ °^(xq) Pq = ^ 

To = - ‘Pq 


dx 


= 1 


dTi 


dy 


= 1 


22 


NACA TN 3340 


Then equation (Cl) hecoiiies 


If |i^ is rewritten in the form 


% - “o * '^n 

where, for gases, 

(C2) 

^i.Q « 1 and Ap.^ << 1 

(C3) 

then equation (Cla) reduces to 

2 

(f|) = ^ ‘‘*0 + PQ 

(Clb) 


P 

where terms in (An^) are infinitesimals of higher order as a result 
of conditions (C3). 

Generally A|j.^ is not constant. Then may be represented by 

the power series 


= 5Z%n'' 

' V =0 

(c^) 

> 

0 

1 

> 

‘WJ 

II 

(C^a) 

where 

, _ 1 /d'’u\ _ + 1 /<a''u\ 

= - ~\7?L 

and 

^^o -^0 

(C5) 

^Tl = fl v" 

V=1 
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Also, let the ray traces "be represented by the power series 

n - r 

cr=0 


where 




l.(^\ s 


(C6) 


(C7) 


Utilizing the assmed series (C4) and (C6) results in the following 
solution of equation (Clb): 


T] = tan Pq ^ sec'^'l 

CT— 1 • • 


00 

sec'^ <P0 

CT— 2^4-^S^ •• • 


(C8) 


where the coefficients c* are functions of the coefficients b^j 
namely 


=t= 1 


ct = i b-: 


* 1 V 

C2 = 2 ^1 


* 1 1, -U 

= 12 ^ 1^2 


and the coefficients c^ are given hy 


Cq = 0 


because n = 0 at ^ = 0, and 


(C9) 


j c* tan Pq sec*^”^ ?>o 
[c* sec (Pq 


(a odd) 
(cr even) 


(CIO) 


Expressions (C9) are obtained by substituting series (C4) and (C6) in 
equation (Clb) and eqmtlng coefficients of like powers of Coeffi- 

cients c* (a even) are identical to the coefficients c^ defineci in 

reference 11 for the case of an axial polnt-light-source. 


Solutions of equation (Clb) for special light-source geometries are 
readily obtainable from equation (C8): 
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(1) The identity (p^q) = (0,0) represents an axial-point light 

soirrce. Then Pq “ "^0 “ equation (C8) reduces to the ray-trace 

equation^, which was derived in reference 11. 

(2) The identity (p,q) = (p^O) represents a narrow slit-light- 
source parallel to the x-axis. Then Pq = 0^ and the resulting siinpli- 

fication of equation (C8) is obvious. 

(3) The identity (p^q) = (0,q) represents a narrow slit-light- 
source parallel to the y-axis. Then Tq = 0^ and the resulting simpli- 
fication of equation (C8) is again obvious. 

It should be noted that equation (C8) increases in complexity in the 
order (l), (2), (3) of the preceding light-so\xrce geometries. 


Optical Distortion 

Because the general fringe-shift equation (3) is a function of the 
optical distortion D, the distortion will be considered first. 

According to equation (Cl) individual ray traces in the field |i(y) 
are refracted in planes. Then^ from figure 2, 


u = 0 


and equation (5) reduces to 


D = H - 

However^ as seen from figures 9(a) and (b) which depict ray-trace pro- 
jections in the xz- and yz-planes, respectively, 

H - D + KL tan Pq = KL tan p^ 

where p = p^ at ^ = L; therefore^ 

D = H - KL (tan P^ - tan Pq) (Cll) 

By definition, H = ^ so that in series form H 

is given by series (C8) with ^ replaced by L, and tan p^^ given 

by the first derivative of series (C8) with ^ replaced by L. Per- 
forming the appropriate substitutions in equation (Cll) and reducing 
the resulting expression, the equation for D, expressed in series form, 
becomes 
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D 


mo 

= ^ 1^0 


(Clla) 


CT=1 


where, by definition. 


%,a = 1 + 


The algebraic sign of D has been chosen to be consistent with the 
coordinate y when D is measiired with respect to the distortion-free 
object point. The resulting sign convention with respect to D is the 
reverse of that used in reference 11. 


When (p,q) = (0,0), equation (Clla) reduces to the distortion 
equation that was derived in reference 11. 


Fringe Shift 


When q = q(y) the integral in equation (3) may be rewritten for 
subsequent simplicity in terms of and the coordinates 5, t), 

Then, by successive substitutions 





TL 


VJ' 


1 + tan2 Xq + 


ll/2 




sec <P/ 


PL 


2 


= sec 


oV^"o 


L + 




where terms in axe infinitesimals of higher order. Because 

l-L^ « 1 and test-section windows are disregarded, the following approx- 
imations can he made regarding especially when <P is small: 
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HaweTCr^ 


sec 

sin 9(-l) 
sec ^^(+l) 
sec ^ (+l) 


J8 


sec <Pq 
sin (pQ 
sec <pQ 
sec <Pl 


sec <P, 


(1 + 

^^0 


tan^ P 

J-i 

sec <Pq 


+ tan^ 


» (1 + Ajig) sec <Pq 

Substituting the appropriate preceding quantities in equation (3) yields 


N = 


sec <Pq 


(^0 - M«.) L + 2 r A|i^d^ ■ A|jijjKL - D sin(pQ cos <Pq 


(C12) 


Equation (C12) can be derived in a more exact form by determining the 
exact expressions for the variovis angles <P as functions of <Pq and 

retaining the multiplicative constants involving and |J.q. 

However, the additional work required is excessive and yields additional 
terms which are appreciable only when <P is impractically large. When 
rewritten in series form, equation (C12) becomes 


sec <P I 


N = 


OP oo 


(p. - n ) L + 2 XI 


kXI Y1 

v=l a=l 


ff+1 


/ oo 

- f ± cr^cr^ ^ ^0 j sin (pQ cos (pQ 

^ CT=1 ’ / 


(C12a) 
a 

obtained by equating coefficients of like powers of ^ in the expression 


where the coefficients c^^^ are functions of the coefficients c^ 
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(v — 1^2;3^...) 


The density ratio p* is obtained by solving ecLuation (C12a) for 
H - utilizing the Dale -Gladstone formula, and dividing both sides 

of the equation by . Th'us , 


* 


P 



X 

L sec «Pq 


v=1 a=l 


0+1 


L<^ 


+ 



V=1 


rT=1 







+ KL 


tan PQjsinf Q cos <Pq 


(C13) 


When (p,q) = (0,0), 
responding equations for 


equations (C12a) and (C13) reduce to the cor- 
N and p* that were derived in reference 11. 


The angle <Pq is a fvmction of angles Pq and ttq. Moreover, Pq 
and Tq are related to the coordinates (p,q) of the light source. 
Specifically, 


tan po = 
tan ro = ±|- 


(C14) 


The conclusion stated in the section Extended or misalined light soiirce 
then results from equations (Clla), (C13), and (C14). 

Equations (Clla) and (C13) yield practical evaluation equations 
for D and p*. However, the coefficients by, c^^, and c^^y must 

first be ejq)ressed as fttnctlons of measurable quantities. The coeffi- 
cients Cg. and are functions of b^. Relations permitting the 

expression of coefficients by as functions of measurable quantities 
are derived in appendix D. 
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APPENDIX D 


RELATION BETWEEN MEASURED AND DISTORTIONLESS FRINGE-SHIFT PROFILES 


ASSOCIATED WITH A ONE -DIMENSIONAL DISTRIBUTION p(y) 


The coefficients up to given hy 

^ 1,1 = '' 1,2 = 0 

*^ 2,1 " ‘^2 ^^ 2,2 "" '^1 

"3,1 = "3 "3,2 = 2ciC2 

"4,1 = "4 "4,2 = "I + 2C]^C3 


(Dl) 


Coefficients c^j are given by equations (CIO) , and coefficients c* are 
given as functions of b^ by expressions (C9). The coefficients b^ 

can be related to measurable space derivatives of the measured fringe- 
shift profile by considering the relation between space derivatives of 
the measured fringe-shift profile and derivatives of the corresponding 
profile that would result if distortion were absent. 


Let 

(1) No(y) denote the fringe-shift profile that would be measured 

if distortion were absent 

(2) N^(y) denote the observed fringe-shift profile 

The y component of separation of the profiles Nq and Np is the 

distortion D. The y-coordinate values associated with a given value 
Nq = Nj) are Jq and y^^, respectively, where 

yp = XO + D 

and D = D(yo). In equations (Clla) and (C12a) all quantities, except 

N, that are variables when associated with a single interferogram are 
defined with respect to initial values jq. However, N is associated 

with values yp. If the fringe-shift profile is distortionless, then 
Np = Nq and is associated with values yp. 
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Consider figure 10, which illustrates hypothetical profiles Mq 
and Nq. In a given fixed interval AN, the incremental slopes of the 
profiles Nq and Nq are, respectively, (AN/Ay)y^ and (AN/Aypi)yQ, 

The ratio of the slopes is 

VAyL Van 

^0 ‘'D 



However^ it results from equation (D2) that 


% 

Ay 


1 


+ 


Ay 


Therefore , 


/AN \ _ 1 /AN\ 


In the limit, as Ay^ ■♦■0, 


dN \ ^ 1 /^\ 

■ 1 + I Wyo 


(D3) 


Relations between higher derivatives of the fringe-shift profiles 
may be determined by differentiating equation (D3) with respect to yp. 

Thus, the n"tli derivative is given by the recurrence formula 


'd”N\ 1 /d”~% \ 

W 1 + g "y Vdyg-7 


(D4) 


Derivatives d^N/dy^ are determined from the measured profile Nq, 
whereas derivatives d’^N/dy^ may be expressed as functions of b^ by 
(differentiating equation (C12a) with the final distortion-contributed 
term omitted (because d’^N/dy^ involves no distortion). By direct sub- 
stitutions equations (D3) and (D4) become, respectively. 
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which include the effect of off-axis source points. The coefficients 

b , c , and c are functions of y and can be differentiated. How- 
V a,v 

ever, derivatives of coefficients b^ are functions of other coeffi- 
cients namely 


^ ^ (v+n); 

dyn “ vJ 4 yVfn vJ ^'+1^ 


(D5) 


as a restilt of definition (C5). Comp'utations involving a finite number 
of terms in equations (D3a) and (D4a) and including derivatives of Np 
up to the derivative must only Include coefficients b^ iq) to bjj^ 

in order to be consistent. Equations (D3a), (D4a), and (D5) then yield 
expressions for by as functions of measurable derivatives d’^w/dyj^. 
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APPENDIX E 


LIMITATIONS OF TWO-TERM EVALUATION EQUATIONS ASSOCIATED 

WITH A ONE -DIMENSIONAL DISTRIBUTION p(y) 

Apparent -ray- trace crossing . - The condition associated with 
apparent-ray-trace crossing can he developed hy considering two ray 
traces which possess initial values Yq and yQ + AyQ, respectively, at 

z = 0, where |yQ + ‘^^yo|^|yo |* Suppose distortions U(yQ) and 
D(yQ + AyQ) are associated with the respective traces. If 

yo + ^yo + ^(yo + ^o) l> l^o + ^^(yo)| 

then "crossing" does not occur hecaiise the larger initial and apparent 
absolute y- values are associated with the same ray trace. However, if 

|yo + ^yo + D(yo + ^yo),|^ |yo + D(yQ)| 

then "crossing" occurs because the trace possessing the larger absolute 
initial value possesses the smaller absoliite apparent -value. When the 
equality is satisfied, the apparent -object points are superimposed. By 
subtracting Jq and D(yQ) from both sides of the preceding equation 

and dividing by AyQ, the following inequality results: 


Ay 


1 


(El) 


where 


AD - D(yQ) - D(yQ + AyQ) 

Ay = yo + % - yo = ^yo 


A practical criterion for avoiding apparent -ray-trace crossing may 
be derived from inequality (El). The limiting condition for crossing 
occurs as Ay 0. Then, inequality (El) becomes 
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The 


derivative of D can be calc\ilated from equation (Clla). Thus, 


^1 

dyl 


f" K 

^ ’^1,0 dy 


which for the case of an axial point-light-source becomes 


® I 

dy| 


t 


dc 


K, 


S2,4,6,... 

Corresponding to the two-term approximation, 

dD 


a 3-0 


(E2) 


dy 


= 0 


because eqmtion (E2) is ultimately independent of b^. However, in the 
three-term approximation, equation (E2) becomes 


1 dD I ^ "K T 2 -L ^ TT ^ 

^ 1 = %, 2 ^ 2 ^ + 2 ^ 1 , 4 ^ 2 ^ 


(E3) 


With convergence assumed, equation (E3) indicates that |dD/dy| = 0 to 
a first approximation when K = l/2. Also, equation (E3) yields a 
criterion for avoiding crossing. If the coefficients b^ (v > 2) are 

assumed to be insignificant, then |dl)/dyl= 1 when 


^1,2'^2^' 


+ i K 


1,4° 2^ 


< 1 


(E4) 


When apparent -ray- trace crossing exists, the experiment must be 
repeated using another value of K beca\ise none of the evaluation 
methods apply. Thus, it becomes desirable to express inequality (E4) 
in terms of quantities that might be estimated before the experiment 
is performed so that a value of K which will avoid crossing can be 
selected initially. The most useful result is obtained by expressing 
b2 as a function of p, namely 


b2 = 


2 ay2 


which results from the Dale -Gladstone formula and definition (CS). Then, 
inequality (E4) immediately becomes the crossing-criterion inequality 
given in the section Apparent -ray-trace crossing . 
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Series convergence . - The series remainders associated with the 
one -term approximation are given by the additional term contained in 
the two-term approximations. The remainders associated with the two- 
term approximations of p* and D are given by the additional terms 
which are contained in the three -term approximation derived in appendix 
F. However^ because in the third approximation the expression for bj_ 

becomes modified by the multiplicative factor additional 

term involving bb^ arises in the remainder equations. The remainders 
associated with the two-term approximation are^ therefore, given by 


where 


|e,(p*)| s 


3kp 




) 






1 ■ 

■ ’^l 


dN 

1,2 = 

<iyD 


dW 

’1,3 = 

dyD 


X 

L 


bg = + — — 


dyg ^^2 


X 

L 


2L 


In general, 

= T - OK 

where t and cr are Integers, and the selected values of the combina- 

p 

tlons b£ ^'t'2 b^ correspond to the largest absolute values 

(denoted by bars) in the interval Yq ^ y = Yq + H traversed by a given 
ray trace. The constants and Xg may be calciilated according to 

the method described in appendix F. 


The present analysis can be applied to develop re si dual -error 
expressions associated with the Wachtell-DeFrate method also. 


34 


NACA TN 3340 


Applicability criterion . - The profiles h2L^ shown in figure 

11 as a function of (cf. appendix F) are approximately straight 

lines with a slope of unity in the vicinity of the origin. Thus, 

• 1, from which it follows that jb 2 j<<l. However, for ^jb2L2j 

\ V’^l/ 

sufficiently large, the profiles deviate drastically from a straight 
line thus indicating that |b 2 j increases considerably. In fact, when 

d^N/dy^ < 0, increasing values of lb L^l eventually correspond to de- 
D 1 2 ' 


D 

creasing values of 




; this condition is absurd from practical 


considerations. The points at which the slopes of the profiles become 
zero (d%/dy2 <o) may be considered, somewhat arbitrarily, as the 
extreme limits for which "true" values of b 2 are obtained. The de- 
fined limits correspond to fixed percentage deviations of the ratio 
x^/x2 from unity, namely 

100^ when K = l/3 

30^ when K = l/2 

60^ when K = 1 

Let the preceding percentage deviations determine the practical limits 
of the profiles for d%/dy| > 0, also. Then because 


^2 V t 2 , d% XL 
dy^ 


the indicated limits correspond to 


-t < ± ^ < 7 when K = ^ 


fl2fj Tlt 

-0.8 < i ^ < 0.8 when 


-0.4 < ± 


dy2 2 

d2w L 
2 2 


K = i 




D 


<0.3 when K = 1 
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The inequalities listed in the section Applicability criterion result 
immediately from the preceding inequalities. 

Test-section windows . - If the test section is bounded by windows, 
then mediums (-1) and (+l) correspond to the initial and final window, 
respectively, and the external atmosphere is contained in mediums (-2) 
and (+2). It has been demonstrated in reference 11 that, associated 
with an axial point -soTirce, the additional optical -path differences con- 
tributed by the presence of windows are: in the final window, 

Pgt [sec <p (+l) - ij (E5) 

and, in the atmosphere succeeding the test section, 

11 ^ [KL - Z(t) + t] [l - sec <p (+2)] (E6) 

where Pg is the refractive index of the window and t is the window 
thickness. The axial separation Z(t) of the apparent- and real-object 
planes is given by 



for small angles , but is exact according to the initial assumption 
(l) in the section GENERALIZED EQUATIONS. 

If off-aocis source points are considered, then contributions (E5) 
and (E6) are replaced by 

Pgt [sec (f (+1) - sec (pj, (+l)J (E8) 

and 

respectively, where Z(t) has been replaced by the right side of expres- 
sion (E7 ) . The initial window introduces no additional contribution to 
the optical -path difference because associated ray traces are parallel 
within the window. 

In reference 11 it has also been shown that terms (E5) and (E6), 
with the quantity KL omitted, cancel to a high order of approximation. 
With terms involving KL disregarded, the only possible additional terms 
of appreciable magnitude contributed by the window are, therefore, 
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Hgt [1 - sec (Pj. (+l)j 


(ElO) 


and 


— t [sec (+2) - 1] 

M-cr 


g 


(Ell) 


These terms, which result from off-axis source points, also cancel if 
sec (+l) and sec (+2) are of the same order of magnitude as 

sec (p (+l) and sec <p (+2). The preceding cancellations result when 
<pQ is zero, in which case the ray traces are given by 

sec 6 = — *■ 

^"0 

If (Pq is not zero, then the ray traces are given by 

1/2 


w]ixch is sinply ecjuation (Cla) rewritten. Thus, the preceding cancel- 
lations result when 




and 


tan^ Tq • 0 


or, in other words, when 



and 



The only remaining term of appreciable magnitude is then 

Pg^KL [sec (pj, (+2) - sec <p (+2)] 


which is actually the last term in equation (3). 
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APPENDIX F 


THREE-TERM EVALUATION EQUATIONS ASSOCIATED WITH A 
ONE -DIMENSIONAL DISTRIBUTION p(y) 

p 

Assuming that = Lq + hj^ri + ^^ 2 ^ ’ evaluation equations 

associated with the three-term approximation are 

(1) P* = 1 + ^ (n J - I Kg^sh^L^ - — Kg^gh^hgL^) 

(2) y = yo = Md " ® 


(3) D = i I ± ^ 


(4) = 


dN 


<iyD 


><2 L 


’'E , ^2 ^ d%XL 


(6) - 1 + + g ^ 


• 2^2 


(7) Xg = 1 + I K 2 ^ 3 b 2 L^ + ^ K2^5(b2L^) 


■2n2 


The evalmtion proceduire is: 

( 1 ) Plot N as a function of y^ graphically from measured data. 

(2) At each datum point measiare slope dW/dyj) of the profile Np. 

(3) Plot dN/dyj) as a function of yj^. 

(4) At each datum point measure slope d^N/dy^ of the profile 
dN/dyp. 

(5) For each datum point con^iute (x2/x^)h2L^. 
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(6) Select appropriate plot of ^2^^ ^ function of (x2/x^)b2L^ 

from figure ll(a) or ll(b), depending upon the algebraic sign of 
d^N/dy^, and determine the value of h which corresponds to the 
appropriate value of K and to the computed value of (x2/x^)b2L^. 

(Plots of (x2/x^)b2L^ as a function of h^L^ are given in figiores 
11(a) and (b) for K = l/3, l/Z, and 1 and for d%/dy| >0 and 
d^N/dyg < 0.) 

(7) For each datimi point conipute b2 = (h2L^)/L^. 

(8) For each dat\jm point compute x^ and X 2* 

(9) For each datum point compute b^, then p*, D, and yQ. 

The density ratio p* as a function of y represents the desired 
profile . 

The measured quantity N and the con5)uted quantities dN/dy^, 

d^N/dyg are associated with the measured profile and values yj). 

All other quantities which are functions of y are associated with 
distortionless values yQ = y. 

The profiles (x2/x^)b2L^ as a function of b2L^ were obtained by 
choosing values of b2L^, computing x^^, X2, and then (x2/x^)b2L^ for 
a given value of K. Thus, profiles can be readily determined for other 
values of K. 
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APPENDIX G 


EQUATIONS FOR CALCULATING PROFILE Np ASSOCIATED WITH A 
HYPOTHETICAL EXPONENTIAL -DENSITY PROFILE 

In order to calculate the profile Nj) associated with an assumed 

hypothetical profile p(y)^ the assumed profile p(y) must he selected 
such that 


_ dll ^ 

pyo+H 

Gy 

,0 , 

.0 ^ 

fn(y) 1 

l2 

- 1 



[^(y) = ^(yo)] 


(Gl) 


and 


OL 




which are involved in computing Np, are readily integrahle in exact, 
or nearly exact, form. Among possible assimied distributions p(y) an 
exponential ftinction, which corresponds to the refractive -index dis- 
tribution 


- (p« - |j^)e"yA 

where p = |j.^ at y = 0^ and a is an arbitrary constant, permits the 

required integrations and corresponds vaguely to a possible boundary- 
layer profile. 

With an exponential density distribution and an axial point source 
assumed, the ray-trace displacements H at ^ = L are given by the 
solution of equation (Gl), namely 

H=alni(l+ cosh t) (G2) 


t 


= 


e 


-yo/2a 


where 
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The ray traces satisfy 

y = yo at ^ = 0 

y = Yl ^ = I' 


where 




At y = Yq, 

v(7o) - I*. - (l*. - (G3a) 

At y = yL, 

^i(yL) = 

The fringe shifts and distortions are given hy equations (C12) and 
(oil), respectively, which become in the present instance 

2 

N =[n. - (1 - K) n(yo) - K liCyp)]^ - 2 ^ tanh(l) (G4) 


D = H - ^2[p(yL) - n(yo)] KL (G5) 

In order to reconpute p*(y) , values of the derivatives d’^N/dy^ 
must be determined. Theoretical values of d’^N/dyg can be computed. 

Derivatives d^N/dy^ are obtained by differentiating eq^uation (G4), 
with q(yo) and dCyp) replaced by the right side of equations (G3a) 

and (G3b), respectively. The required first two derivatives are given by 



(G6b) 
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where 



Derivatives d’^H/dy^ follow from eqToation (G2). Thus, 


where 


sinh t ^ 

M = a ^ 

dy 1 + cosh t 



(S) 


2 d^t 

+ sinh t — ?r 


a 


1 + cosh t 




d% 

dy^ 


2a 


u 




D 


sech^ % + 


dt -t 
dy 2a 

d^t _ 
dy^ 4a^ 


(G7a) 

(G7h) 

(G8a) 

(G8h) 

(G9a) 

(G9h) 

(G9c) 

(GlOa) 


(GlOb) 
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Derivatives d’^N/dyg are functions of derivatives d^N/dy^ given 'by- 

equation (D3) and the recurrence formula (D4) . The required first two 
derivatives of D are obtained by differentiating equation (G5) follow- 
ing the previo-usly mentioned substitutions for l-i(yo) l-'’(yL)* Then, 


dy 


dy 


^ ^ + 2aK 


d% 

dy2 


(Glia) 


^ ^ + 2aK ^ (Glib) 

dy^ dy2 dy^ 

The profile Np is determined by computing equations (Gl) to (G5), 
and the derivatives d°N/dyg (n = 1,2) are calculated from equations 

(G6) to (Gll). 


The selected values of the constants were 


= 1.0000792 (p. = 6.77X10"^ slug/ft^) 

u = 1.0000475 (p = 4.06X10"^ slug/ft^) 
w w 

a = 0.0058 in. 

X = 2.15X10“^ in. (mercury green line) 
X = 0.117 ft^sl^ 


which, except for a, are 
example of reference 11. 


identical to the values Involved in the 
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APPENDIX H 


EXPERIMENTAL VERIFICATION OF REFRACTION EFFECT 

The conduction of heat by air between two horizontal surfaces pro- 
vides a readily produced and reasonably predictable density distribution 
which may be utilized for verifying the refraction effect and two-term 
evaluation equations. Consider two identical rectangular-shaped blocks 
placed one a,bove the other and separated by a small distance, with fac- 
ing surfaces A and B corresponding to the bottom surface of the top 
block and the top surface of the bottom block, respectively. Suppose 
that both blocks contain a reference hole which extends throughout the 
horizontal length of the block and is oriented parallel to the optical 
axis of the interferometer. 

When the upper block is heated, a constant temperature gradient is 
produced in the airspace between the blocks, whereas the temperatures 
within the reference holes are constants. If the blocks are good heat 
conductors, then the temperature T^_ within the reference hole through 

the upper block and the temperature T^ at surface A are equal. How- 
ever, (dp/dy)^ = 0 whereas (dp/dy)^ = constant ^ 0, so that light which 

traverses the space between the blocks is refracted whereas light which 
traverses the reference hole is not. Thus, the optical -path difference 
between light passing adjacent to A and light traversing the reference 
hole should theoretically be given by 


These equations contain only contributions attributable to refraction. 

In order to properly associate the refraction and reference fringes, 
the algebraic sign of the fringe shift corresponding to the reference 
region in the lower block must be determined. Because T 2 where 

T 2 is the temperatuire of the lower block, |i 2 ^ so that N 2 ^ 0. 
Therefore, it follows from equation (Hi) that: 



(HI) 


and the distortion should be given by 



(H2) 


44 


NACA TN 3340 


(1) If Yi <Z/ "5, then > 0 and the fringes between the blocks 

are displaced toward the zero-order fringe in the lower reference region 
(fig. 12(a)). 

(2) If K = 2/3, then = 0, that is, no fringe shift occurs 
(fig. 12(b)). 

(3) If K >2/3, then Na < 0 and the fringes are displaced away 
from the zero-order fringe in the lower reference region (fig. 12(c)). 

A corresponding interferogram utilizing unfiltered light will permit 
the zero-order fringe to be distinguished. 

The model shown in figure 13 was constructed in order to simulate 
the aforementioned conditions. The model consists -of two Duraluminum 
blocks 10-inches long which are separated by two 0.14-inch thick, 10- 
inch long insulating strips of glass. The top surface of the upper 
block is radiant -heated by a Nichrome-wire heater consisting of Nichrome 
resistance wire wound with varying pitch around a thin sheet of mica. 

The heating element and most of the upper block are insvilated from the 
external atmosphere by a Trans ite lid in order that the heat -flow rate 
from the upper block to the lower block will be a maximum for any given 
heating current. A reference hole was machined the length of each 
block, and a reference pin was projected into each reference hole for 
focussing purposes and as a measuring reference. 

Figures 14(a) and (b) depict typical interferograms of the heated 
air for K = 1/2 and 1, respectively; these interferograms were ob- 
tained by utilizing unfiltered light from a magnesium spark in conjunc- 
tion with a horizontal slit. It is apparent that the temperature gra- 
dient between the blocks is nearly constant. Curvature of the fringes 
near the upper s\jrface of the lower block is caused by reflection. 

In figures 15(a) and (b), experimentally determined values of 
optical-path differences and distortion, respectively, are compared 
with theoretical values indicated by equations (Hi) and (H2). The 
qmntities N^X and D^ are plotted with respect to the generalized 

coordinate (dN/dyj))j^X which yields a single curve for all wave lengths 
X. Actually, the experimental data include measiirements at two wave 
lengths, namely 4481 A (magnesiimi spark) and 5461 A (mercury arc). 
Agreement of the experimentally determined points with the predicted 
c-urves is regarded as good considering the errors involved in the ex- 
periment and the fact that the difference in values amounts to the 
error of measuring a small error . 

The most important experimental errors probably were end effects, 
thermal expansion of the apparatus, nonuniform heating of the apparatus, 
inclination of reference fringes, and poor definition of the apparent 
location of surface A. 


3197 


3197 


NACA TN 3340 


45 


In attempting to minimize end effects, the model was initially 
hounded hy • 0.8-inch-thick glass windows separated from the model ends 
hy thin cork spacers. The windows were slid into position immediately 
preceding each photographic exposure. The experimental data were then 
in poor agreement with theory. Apparently, insertion of the windows 
prior to the photographic exposure, ev^en for very short periods of time, 
allowed appreciable nonuniform heating of the glass. Air disturbances 
were not evident. It would a’ppear that extreme caution should be ex- 
ercised in interpreting interferograms of heated models when the model 
contacts bounding glass windows. In order to eliminate the window 
effects, the windows were removed, and the ends were left unbounded. 

A new set of measurements rendered the results presented in figures 
15(a) and (b). The interference pattern obtained across the end of the 
model with the open-end arrangement was examined visually while the 
model was heated. The endwise boundary layer was found to be of nearly 
constant thickness (about l/4 inch) in the region of interest, and the 
temperature gradient was essentially parallel to the model length at 
all pertinent temperatures. Thus, the end effects were regarded as in- 
appreciable when the windows were removed. 

Thermal expansion of the model was considered in the vertical di- 
rection, but expansion in the lengthwise direction was neglected. 

Teirperatures were measured throughout the model by means of thermo- 
couples located at the points indicated in figure 13(b). A typical 
temperature distribution is also indicated in the figure, where the 
temperature difference at each point is noted with respect to the upper- 
most thermocouple location in the midspan plane. Inclination of refer- 
ence fringes with respect to surface A, caused by initial misorientation 
of the fringes and the small vertical temperature gradient within the 
upper block, was considered in the measurements. 

The existence of refracted ray traces in the vicinity of surface 
A tended to diminish the definition of the apparent surface location. 
The difficulty of locating surface A, especially when K = l/2, may 
account for the apparent systematic disagreement between theory and 
e^eriment when K = l/2. 
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Figure 1. - Zehnder-Mach interferometer and characteristic ray traces. 
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Figure 2. - Light-path detail within test section. 
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(a) 



(t) 


Figure 3. - Examples of apparent -ray-trace crossing. 
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Figure 4. - Light-path detail preceding collimating lens. 
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Density ratio 



Figure 5. - Re-evaluation of exponential density profile where L = 1.8 inches. 
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(Id) Two - term approximation. 

Figure 5. - Continued. Re-evaluation of exponential density profile where L = 1.8 inches. 


cn 

CM 


MCA TN 3340 


Density ratio 


Oi 



Figure 5. - Concluded. Re-evaluation of exponential density profile where L = 1.8 inches. 
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Figure 6. - Re-evaluatlon of exponential density profile where L = 3.6 inches. 
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(t) Two-term approximation. 

Figure 6. - Continued. Re-evaluation of exponential density profile where L = 3.6 inches. 
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(c) Three-term approximation. 

Figure 6. - Concluded. Re-evaluation of exponential density profile where L = 3.6 inches. 
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Figure 7. - Comparison of re-evaluated density profiles. 
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(b) L = 3.6 inches. 

Figure 7. - Concluded. Comparison of re-evaluated density profiles. 
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Distortion gradient 



Figure 8. - Distortion and distortion gradient associated with object planes corresponding to K = l/s, l/2, 
and 1 and assumed density profile, where L = 3.6 Inches. 
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Figure 9. - Ray-trace projections In presence of a unidirectional 
distribution p = p(y) • 
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Figure 10. - Relation of observed profile Np(y) 
to undlstorted profile NQ(y) . 
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Figure 11 


64 


MCA TN 3340 





Figure 12. - Fringe shifts associated with various focal planes in hot-plate model 
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(b) Schematic diagram. (Temperature drop indicated with respect to 
temperature ( 272 ^) at uppermost midspan thermocouple location.) 


Figure 13. - Concluded. Hot-plate model. 
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Concluded. Interferogram of heated model 
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Figure 15. - Verification of refraction- error formula associated with 
constant density gradient. 
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(b) Distortion; K = 1. 

Figure 15. - Concluded. Verification of refraction- error formula associated 
vith constant density gradient. 
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